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Abstract 

The two known non-perturbative theories of localization in disordered wires, 
the Fokker-Planck approach due to Dorokhov, Mello, Pereyra, and Kumar, 
and the field-theoretic approach due to Efetov and Larkin, are shown to be 
equivalent for all symmetry classes. The equivalence had been questioned as a 
result of field-theoretic calculations of the average conductance by Zirnbauer 
[PRL 69, 1584 (1992)], which disagreed with the Fokker-Planck approach in 
the symplectic symmetry class. We resolve this controversy by pointing to an 
incorrect implementation of Kramers degeneracy in these calculations, and we 
derive modified expressions for the first two conductance moments which agree 
well with existing numerical simulations from the metallic into the localized 
regime. 

PACS numbers: 72.15.Rn, 72.10.Bg, 02.50.-r 



1 



Typeset using REVTeX 



I. INTRODUCTION 



There are two known approaches to the theory of phase- coherent conduction and local- 
ization in disordered wires: The first is the Fokker-Planck approach of Dorokhov, Mello, 
Pereyra, and Kumar. The second is the field-theoretic approach of Efetov and Larkin, 
which leads to a supersymmetric nonlinear a model. ii Both approaches provide a descrip- 
tion of quantum transport that is independent of microscopic details of the disordered wire. 
The only properties which enter are its length L, the elastic mean free path £, the num- 
ber N of propagating transverse modes at the Fermi level (referred to as "channels"), and 
the symmetry index (3 G {1,2,4} (depending on the presence or absence of time-reversal 
and/or spin- rotational symmetry). In the first approach, the transfer matrix is expressed 
as a product of a large number of random matrices. As more matrices are added to this 
product, the transmission eigenvalues T n execute a Brownian motion. (The T n are the N 
eigenvalues of the transmission matrix product tH.) The resulting Fokker-Planck equation 
for the L-dependence of the distribution P(T 1; . . . , T/v) is known as the Dorokhov-Mello- 
Pereyra-Kumar (DMPK) equation. In the second approach, one starts from the random 
Hamiltonian of the disordered wire and then expresses averages of Green's functions!'! or 
moments of the transmission eigenvaluesi^ as integrals over matrices Q containing both 
commuting and anticommuting variables. These so-called supermatrices are restricted by 
the nonlinear constraint Q 2 = 1 and give rise to a field theory known as the one-dimensional 
nonlinear a model. 

In the last decade, research on the Fokker-Planck and field-theoretic approach has pro- 
ceeded quite independently. Recently, exact results for the average conductance (G), its 
variance varG, and the density p{T) = (I] n <5(T — T n )) of transmission eigenvalues were 
obtained from both approaches. For the unitary symmetry class (no time-reversal sym- 
metry; (3 = 2), the DMPK equation was solved exactly by Beenakker and Rejaei. i The 
construction of a set of biorthogonal polynomials for this exact solution then allowed for 
the exact computation of (G), varG, and p(T) for arbitrary N and L in the case (3 = 2.0 
Although there exists a formal solution for the other two symmetry classes [orthogonal class 
(time-reversal symmetry without spin-orbit scattering; (3 = 1) and symplectic class (time- 
reversal symmetry with spin-orbit scattering; (3 = 4)],EJ no exact results for (G), varG, and 
p(T) have been obtained. Concerning the a model, an important and substantial progress 
was the development of "super Fourier analysis" by Zirnbauer.ll2l This allowed the exact 
calculation00 of (G) and var G for all (3 in the thick-wire limit N — > oo, L/i — > oo at fixed 
ratio Nl/L. The eigenvalue density p(T) was computed from the a model by Rejaei,0 in 
the thick- wire limit and for the case (3 = 2. 

If one takes the thick- wire limit of the (3 = 2 results for (G), varG, and p(T) from the 
DMPK equation, they agree precisely with those from the a model.0'0 For (3 = 1 and 4, a 
comparison of the two approaches has only been possible in the metallic regime f « L « NL 
where the results for (G) and var G from the DMPK equation! - ! and from the a model!0Jll 
agree with conventional diagrammatic perturbation theory.llUHl The equivalence of the two 
approaches outside the perturbative regime has been questionedlll as a result of recent work 
by Zirnbauer,0 and by Mirlin, Muller-Groeling, and Zirnbauer.0 Starting from the a model 
in the thick- wire limit, they obtained a finite limit (G) —>■ e 2 /2h as L/N£ — > oo in the 
case (3 = 4. On the other hand, one can prove rigorously^! that the DMPK equation gives 
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Hindoo (G) = for all (3. It was this puzzling contradiction which motivated us to search 
for a general proof of equivalence of the DMPK equation and the a model, without the 
restriction to (3 = 2. 

In this paper, we present a general proof of the equivalence of the two approaches, which 
applies to all three symmetry classes (3, to all length scales L, and to the complete dis- 
tribution of transmission eigenvalues described by the p-point functions p p (Ti, . . . , T p ) = 
(N\/(N — p)\) J dT p+ i . . . J dTjsi P(Ti, . . . , T/v) for arbitrary p. We cannot relax the assump- 
tion that the number N of propagating channels in the disordered wire is ^> 1, since it 
is needed for the derivation of the one-dimensional a model.0 However, we can consider 
the a model formulation of a thick disordered wire which is coupled to the leads by means 
of a point contact with N\ < N transmitted modesjl and show that it is mathematically 
equivalent to a DMPK equation for a wire with N% propagating channels. The equivalence 
proof demonstrates that lim^^ 00 (G') = in the a model, in apparent contradiction with 
Zirnbauer's work. We have reexamined the calculation of Refs. |H] and [TT], and argue that 
for (3 = 4 the Kramers degeneracy of the transmission eigenvalues was not taken into account 
properly in the super Fourier analysis. This leads to a spurious "zero-mode", which does 
not decay as L — > oo. Restoring Kramers degeneracy, we obtain modified expressions for 
(G) and var G which decrease exponentially in the localized regime and moreover agree well 
with numerical simulations.0 

Both the a model and the DMPK equation were derived from a number of different 
models for a disordered wire. The original derivation of the DMPK equation by Dorokhov,0 
which started from a model of N coupled chains with defects, was followed by the random- 
matrix formulation of Mello, Pereyra, and Kumar.! These authors considered a product 
of random transfer matrices, drawn from an ensemble of maximum entropy. Later it was 
shown that the DMPK equation is insensitive to the choice of the ensemble, the only relevant 
assumptions being weak scattering (mean free path I much greater than the Fermi wave 
length \p) and equivalence of the scattering channels.011 It is this latter assumption which 
restricts the DMPK equation to a wire geometry. From the mathematical point of view, 
the DMPK equation is the diffusion equation on a certain coset-space of transfer matrices£j 
The one-dimensional a model was originally derived by Efetov and LarkirJH from a white 
noise model for the disorder potential. Two later derivations used random-matrix models 
for the Hamiltonian of the disordered wire. Iida, Weidenmuller, and Zuk (IWZ) adapted 
Wegner's n-orbital modelil to the study of transport properties.i In this description, the 
wire is modeled by a large number of disordered segments in series, each segment having a 
random Hamiltonian drawn from the Gaussian ensemble. An alternative derivation of the a 
model, due to Fyodorov and Mirlinjli uses a random band matrix to model the Hamiltonian 
of the disordered wire. In the present paper we follow Ref. [11] and use the IWZ formulation 
of the a model. 

Our proof of equivalence of the DMPK equation and the a model builds on the ideas 
which were used by Rejaeill! to calculate p(T) from the a model for (3 = 2. Inspired by 
Nazarov's diagrammatic calculation of p(T) in the metallic regime,@ Rejaei introduced a 
generating function F which depends both on the transmission eigenvalues T n and on the 
radial parameters 9i of the supermatrices in the unitary a model. Rejaei was able to solve 
the Id a model exactly for (3 = 2 and thus obtained the density p(T) as a function of L by 
taking derivatives of F with respect to the 0j's. The resulting p(T) could then be compared 
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with the result from the DPMK equation!! We introduce a more general generating function 
which allows us to establish the equivalence of the a model and the DMPK equation at the 
level of p-point functions p p (Ti, . . . ,T p ), without actually having to compute this function. 
This approach works also for (3 = 1 and 4, where no explicit solution of the a model is 
available. 

The outline of the paper is as follows: In Sec. 0, an outline of the equivalence proof 
is given. The full proof for the a model with 8x8 supermatrices follows in Sees. [TIT] and 
IV| , with technical material in Apps. 0-0. For the p-point functions p p (Ti, . . . ,T p ), we 
have to consider the a model with 8p x 8p supermatrices. This extension is described in 
App. [D]. In section [V], we discuss the symplectic symmetry class (/3 = 4) in relation to Refs. 
[T~0| and |ll[ By accounting for Kramers degeneracy, we obtain modified expressions for (G) 
and varG, which we compare with numerical simulations of the IWZ model by Mirlin and 
Muller-Groeling.lH We conclude in Sec. |V]. 



II. OUTLINE OF THE EQUIVALENCE PROOF 



Although our equivalence proof is technically rather involved, the basic idea can be 
described in a few paragraphs. In this section, we present an outline of the equivalence 
proof for the small a model (8x8 supermatrices). The details are given in the following two 
sections and in the appendices [A] — |C[ Appendix [D| contains the necessary modifications to 
extend the proof to a models with supermatrices of arbitrary size. 

Part of the complexity of the problem is that the a model and the DMPK equation focus 
on totally different objects. In the a model, transport properties are expressed as functional 
integrals over supermatrices (A supermatrix is a matrix containing an equal number 

of commuting and anticommuting elements. We follow the notation and conventions of Refs. 
[8], H and [11].) For the small a model the 8x8 supermatrices are parameterized asBi 



Q = T _1 AT, A 
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where u and v are pseudo-unitary 4x4 supermatrices. Notice that Q satisfies the non-linear 
constraint Q 2 — 1, hence the name "non-linear" a model. (The letter a is used for historical 
reasons.) The 4x4 supermatrix 9 is called the radial part of Q. It has the form 
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While the a model works with the radial part of a supermatrix, the DMPK equation 
works with the radial part of an ordinary matrix (containing only commuting elements). 
This is the transfer matrix X. The radial part of X is an N x N diagonal matrix A, related 
to the eigenvalues of XX'. The eigenvalues of XX' come in N inverse pairs e ±Xn , related 
to the diagonal elements A n of A by A n = sinh 2 x n . For (3 — 4 the eigenvalues are twofold 
degenerate (Kramers degeneracy). The matrix A then contains only the N independent 
eigenvalues. The conductance G is directly related to the A n 's by the Landauer formulai@ 



G 



2^ v 
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n=l 



2e^ 
h 



N 



n=l L 



(2.2) 



since the N independent transmission eigenvalues T n are related to the A n 's by T n = (1 + 

We connect both approaches by considering a generating function F(9, A) which depends 
on both radial matrices: 



iv 



F{0,X) = l[f{6,\ n l 

n=l 

f{9, A) = Sdet^ /2 [A + cosh 2 (#/2)" = 

(1 + 2A + cos(fl 3 + 9a)) (1 + 2A + cos(g 3 - gj) " 2 
(1 + 2A + cosh(^i + 2 )) (1 + 2A + cosh(^i - 9 2 )) 

d=lif/3 = l,2; d = 2 if /? = 4. 



(2.3a) 
(2.3b) 



(2.3c) 



The symbol Sdet stands for the superdeterminant of a supermatrix. For (3 = 2 this is the 
generating function introduced by Rejaei. 

An ensemble of disordered wires of length L provides a distribution of A. The ensemble 
average (F(9, A)) contains all statistical properties that are accessible from the small a 
model. These include the average conductance (G), its variance varG and the density of 
transmission eigenvalues p(T). We explain in appendix [A] how to extract these quantities by 
taking derivatives of (F(9, A)). The average (F(9, A)) can be determined by each of the two 
approaches independently, in terms of a partial differential equation for the L-dependence 
and an initial condition at L = 0. For the a model on the one hand, the evolution equation 
reads 



F(9, A) 



(2.4a) 



where is the (radial part of the) Laplacian on the a model space, and where £ = /3N£ is 
the localization length. The explicit form of Ag is given byQ 



(2.4b) 



where the sum runs over the independent coordinates 9i [see Eq. (|2.1d )] and J{9) is the 
integration measure for the radial decomposition (|2.1|), 
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J{Oi, 02, 9 3 ) = sinh 9 1 sinh 9 2 sin 3 6 3 II sinh " 2 (l^i + s ^ + is ^)) if /? = 1, 

Sl,S2=±l 

J(e l ,9 3 ) = sinh 6 1 sin # 3 J] sinh -2 (±(0i + isi0 3 )), if /3 = 2, (2.4c) 

31=±1 

J(e 1 ,e 3 ,9 4 ) = sinfl 3 sinfl 4 sinh 3 ^ J] sinh" 2 {\{6 X + zs x 3 + is 2 # 4 )) if /? = 4. 

Sl,S 2 = ±l 

The DMPK equation on the other hand, yields the evolution equation 

JL(f0,\)) = ^(d x F0,\)), (2.5a) 

where D~ x is a second order differential operator in the parameters A n , 

- N d d 
Di = J~\\) E or ^(A) An(l + A n ) — , (2.5b) 

n=l n n 

J(X) = II |An-A m |^. (2.5c) 

The key ingredient of the equivalence proof is the identity 

A $ F0,\)=D X F0,\), (2-6) 

which shows that the evolution with L of (F0, A)) is the same in both approaches. Showing 
that the initial conditions at L = coincide as well, completes the equivalence proof. 

III. ONE-DIMENSIONAL a MODEL 

We begin the detailed exposition of the equivalence proof with a formulation of the a 
model. As in Ref. 11, we use the formulation of Iida-Weidenmiiller-Zuk (IWZ).i 



A. The IWZ model 

The IWZ mo deffl applies Wegner's n-orbital mo deli to a wire geometry and sup- 
plements it by a coupling to ideal (not disordered) leads, as in Landauer's approach to 
conduction.il The left and right leads (labeled by indices 1 and 2) contain Aq and N 2 prop- 
agating modes each (per spin direction for (3 = 1,2, or per Kramers doublet for (3 = 4). 
The disordered wire of length L is assumed to consist of K segments in series (Fig. |l|). The 
Hamiltonian H of the disordered wire without leads is represented by a matrix H^ u , where 
the upper indices i, j label the segments 1 < i, j < K and the lower indices /i, v label the M 
states (per spin direction or Kramers doublet) within each segment. The elements of H are 
real ((3 = 1), complex ((3 = 2) or quaternion ((3 = 4) numbers. The coupling between the 
states inside one segment is described by the matrices H", which are distributed according 
to the Gaussian ensemble 

P(H li ) = const, x exp (-^Mv^Tr (H 11 ) 2 ) . (3.1) 
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Here v\ is a parameter which governs the level density at the Fermi level (E = 0). The cou- 
pling between the states of adjacent segments is given by another set of Gaussian distributed 
random matrices = (fP 1 )' (with coupling parameter v 2 ), 

P(H ij ) = const, x exp (-\(3M 2 v^ 2 Ti H ij H ji ) , 

j = i±l. (3.2) 

Segments which are not adjacent are uncoupled, — if \i — j\ > 2. The coupling to 
the ideal leads is described by a fixed KM x (Ni + N 2 ) rectangular matrix W = W\ + W 2 
with real {(3 = 1), complex (/3 = 2) or quaternion (/3 = 4) elements. The matrix W has 
elements W l , where i labels the segment, \i the states in the segment, and n the modes in 
the leads. The elements of W\ (which describes the coupling to lead 1) are nonzero only for 
i = 1 and 1 < n < Ni\ the elements of W% (coupling to lead 2) are nonzero only for i = K 
and N 1 <n<N 1 + N 2 . 

The scattering matrix S (matrix elements S nm ) of the system at energy E is given byi 

S = 1 -2iriW^(E - H + ittWW^)- 1 W. (3.3) 

The indices n, m correspond to lead 1 if 1 < n,m < N± and to lead 2 if Ni < n,m < Ni + N 2 . 
The reflection and transmission matrices r,r',t,t' are sub matrices of S, 



S = 1 ,) ■ (3-4) 




Since S is unitary, the products tH and t'H' have the same set of non-zero eigenvalues, 
denoted by T n = (1 + A n ) _1 . (If N 2 > N\ there are also N 2 — Ni transmission eigenvalues 
which are zero, and can therefore be disregarded.) 



B. The generating function 

We now define the generating function E(6, A) introduced in the previous section. We 
start from the the relationship ( p.3| ) between the scattering matrix and the Hamiltonian in 
the IWZ model. We consider the generating function 

F = SdeH(£ _ n + nxW x W\Q + iwW 2 W^A), (3.5a) 
H = Hl s if P = 1,4; H= (ReH)l 8 + i(lmH)r 3 if (3 = 2. (3.5b) 

Here 1§ is the 8x8 supersymmetric unit matrix and r 3 is a diagonal matrix with elements 
(1, — 1, 1, — 1, 1, — 1, 1, — 1). The matrix A was defined in Eq. ( 2.1a| ). Note that Q is an 



arbitrary supermatrix as in Eq. ( |2.1|) and that it replaces the matrix A in the coupling term 
of lead 1. In App. [TJ we show that F depends only on the radial part 9 of the matrix Q and 
that the only dependence on H is through the transmission eigenvalues T n = (1 + A n ) _1 . We 
also show that Eq. ( |3.5| ) reduces to the function F(9, A) defined in Eq. ( |2.3| ) of the previous 
section. 
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In the following, we evaluate the ensemble average (F) using the supersymmetric for- 
malism. We first express (F) as a Gaussian integral over an 8Mi^-dimensional supervector 

(F) = Q Vi) exp (§z ^A(E -H + ntW x W\Q + mW 2 W\k + teA)ip) \ . (3.6) 

The convergence of the Gaussian integral is assured by the parameterization ( |2.1| ) of the 
matrix Q. Performing the standard steps, described in Refs. ^ and [11, we obtain in the 
relevant limit M — > oo 

(F) = JdQiJ dQ K /i(Q,Qi) / 2 (A,Q x )W(Qi,Qjf), (3.7a) 
Qk)= JdQ 2 ...J dQx-i exp ^-^f ]T Str(Q, Q i+1 ), ) (3.7b) 

=exp f-^EStr ln(l + x„gg 1 )j , (3.7c) 

/ i N1+N2 \ 
f2{Q,QK)=exp[--d ]T Strln(l + a; n QQ^) . (3.7d) 

\ 2 n=Ni+l J 

The numbers x ra denote the eigenvalues of the matrices (ir/v \)Wl W\ (if 1 < n < iVi) or 
(Tr/uOWjWa (if iVi < n < iVi + iV 2 ). The integer d was defined in Eq. (£3cj). 

Following Ref. |ll|, we consider the limit v\ <C u|. Then the sum in ( [3.7b] ) can be replaced 
by an integral and the (Jj-integrals yield a path integral. The discrete number of segments K 
becomes the continuous (dimensionless) variable s. The propagator ( |3.7b|) can be identified 
with the heat kernel of the supersymmetric space, determined by the heat equation0 

2P(v 2 /v 1 ) 2 -?-W(Q\Q") = A Q ,W(Q',Q"), 

\imW(Q',Q")=6(Q',Q"). (3.8) 

The precise definition of the Laplacian Aq and the detailed justification of Eq. (|3.8| ) are 
contained in Ref. [11] (Aq in Eq. ( |3.8| ) differs by an additional factor @/(8d) with respect to 
the notations of Ref. ]TT| ). We thus arrive at the expression 

(F) = JdQ'J dQ" h{Q,Q')W{Q',Q")f 2 {K,Q"). (3.9) 

The next step is to notice that fi{Q,Q') has the same symmetry as the heat kernel, 
i.e. /i(T _1 QT, T _1 Q'T) = fi(Q,Q') where T is an arbitrary element as described in fl2.1a| ). 
This implies Aq//i(Q, Q') = Aq/i(Q, Q') and hence (F) also satisfies the heat equation 

2P(v 2 / Vl ) 2 j- s (F) = A Q (F). (3.10) 

Since (F) only depends on the radial part 9 of Q, it is sufficient to consider the radial part 
Aq of the Laplacian Aq. This radial part A^ can be written as in Eq. ( |2.4b| ). We thus find 
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that the ensemble average {F(9, A)) of the generating function defined in Eq. fl2.3|) satisfies 
the partial differential equation 

2(3(v 2 /v 1 ) 2 -^(F(§, A)) = A S (F(0, A)), (3.11) 



with the initial condition implied by Eq. ( |3.8|) , 

Km(F{§,\)) = J dQ'f\{Q,Q')f 2 {\,Q'). (3.12) 



Together, Eqs. ( |3.11| ) and ( |3.12| ) determine the ensemble average of the generating function 
F(6, A) evaluated in the framework of the nonlinear a model. 

The two limits of the IWZ model which were needed for the derivation of Eq. ( |3.12 ) 



M — > oo and v\jv\ — > 0, restrict the validitv of Eq. (|3.12j) to the case of weak disorder 
(l ^> Xf) and thick wires (N ^> 1) respectivelyBS Whereas the requirement of weak disorder 
is also needed for the DMPK equation, the requirement that the number of channels in the 
disordered wire be large is not. To see how the latter requirement follows from the condition 
v\ <C u|, we consider the expression for the average conductance (G) in the diffusive metallic 
regime (l<L< N£)jM 

2p 2 NP 2p 2 Av 2 
(G) = ^ = (3.13) 
n L h vfs 

Taking the linear dimension of a segment of the disordered wire in the IWZ model of order 
P (i.e. s ~ L/P, see Ref. ^), we find that v\ <C v\ corresponds to > 1. However, no 
restriction has been put to the numbers Ni and N 2 of propagating channels in the leads in 
the above derivation of the a model, which allows us to consider finite values of N\ and A^2- 
This situation corresponds to the case in which the thick disordered wire is coupled to the 
leads 1 and 2 by means of point contacts, with N\, N 2 open channels. As in Ref. |TT], the 
case of a disordered wire without point contacts is recovered in the limit N\, N 2 — > oo. 

We conclude this section with some remarks about the choice of initial conditions. In 
usual a model calculations,0'0'0 one considers ideal coupling (x n = 1, n = 1, . . . , N\ + N 2 ) 
and identifies N = N\ = N 2 (equal number of channels in the leads and in the wire). In the 
thick-wire limit N —> oo the function fi(Q, Q') is just the delta functionS 5(Q, Q'), and (F) 
becomes identical to the heat kernel itself [cf. Eq. (|3.9|) 1: 

(F) = W{Q,A), N 1 = N 2 = N^>1. (3.14) 

For P = 2, this result was derived by RejaeiB In this case (F) has the delta-function initial 
condition lim s ^ {F) = S(Q,A). To make contact with the DMPK equation, we need a 
different "ballistic" initial condition, such that all T n 's are unity in the limit of zero wire 
length. To achieve this, we take ideal coupling and assume that one of the leads has many 
more channels than the other. To be specific, we fix N\ and take the limit N 2 —>■ oo. One 
then finds the initial condition 



lim(F) = exp (-fA^Str ln(l + QA)) 

/ COS #3 + COS U4 \ 



cosh 6>x + cosh 9 2 



1 < JVi < N 2 . (3.15) 
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In the next section, we will see that this is precisely the ballistic initial condition of the 
DMPK equation. 



IV. DMPK EQUATION 

Let us now evaluate the ensemble average of the generating function Q2.3|) from the 
DMPK equation. The DMPK equation is a Fokker-Planck-type equation for the L-evolution 
of the probability distribution P(X) of the A„'s:§Hl 

N a a 

? W n Ul + ^JWg^J-'WPW, (4-la) 

■JW= II |An-A m |^, (4.1b) 

n>m 

where i denotes the mean free path in the disordered wire and N the number of propagating 

modes. There is no restriction to iV ^> 1 in the DMPK approach. We take the ballistic 
initial condition 

N 

limP(A) = Tp(A n -0 + ). (4.2) 

L — >U 

n=l 

The DMPK equation implies for F(9, A) the evolution equation!'! 



— (F(9, \)) = —Jd\ l ...J d\ N F(6, A) P(A) 



= - i ((3N + 2-P)- 1 (D x F(§,\)), (4.3) 

with the differential operator given by Eq. (|2.5b| ). In Appendix y we prove the algebraic 
identity between the two different types of Laplacians ( |2.4b| ) and (|2.5b| ) applied to F(9, A), 

A 9 F(O,\) = D X F0,\). (4.4) 

From Eqs. ( |4.3|) , and ( |4.4|) we conclude that that the average {F(9, \)), calculated in the 
framework of the DMPK equation, also fulfills the evolution equation ( |3.11| ) of the nonlinear 
a model, provided we identify [cf. Eq. (|3.13|) 1 

4 Nf P 

s L pL 

Here we introduced the localization length £ = /3N£ (notice that the definition of £ in Ref. 



lT| differs by a factor 2/ j3). 

It remains to compare the initial conditions. The ballistic initial condition for the DMPK 
equation implies 
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Nd 

I , (4-6) 



which equals the initial condition ( |3.15 ) for the nonlinear a model (The thick-wire limit 



lim i ^ (F) = 5(Q,A) is obtained by letting N — > oo in the above expression). This proves 
the equivalence of both approaches, as far as the generating function (|2.3|) is concerned. In 
Appendix we extend the equivalence proof to p-point functions p p (Ti, . . . , T p ) for arbitrary 
p. 



V. THE CONTROVERSIAL SYMPLECTIC ENSEMBLE 

The main motivation of this work was to resolve a controversy between the DMPK 
equation and the one-dimensional a model in the symplectic symmetry class {(3 = 4). On 
the one hand, the DMPK equation impliesll! (G) —>■ as L — > oo. On the other hand, 
ZirnbauerS finds from the a model that (G) — > |e 2 /h as L — > oo. 

The equivalence proof presented in this paper has as a logical consequence that (G) — > 
as L — > oo if (G) is evaluated in the framework of the a model. To demonstrate this, we 
apply the argument of Ref. [13]. The DMPK equation implies for the average dimensionless 
conductance g = X) n (l + ^n) -1 the evolution equation! 

£^ = -/% 2 >-(2-/3)<<72>, (5.1) 



with g 2 = S n (l + A n ) 2 . This relation also follows from the evolution equation ( |3.11| ) of 



the a model (expanding the generating function for small 9{ and applying the results of 
appendix 0). Since < g<i < g 2 , we have 

^ " -\^ 2) ~ °- (5 ' 2) 

We suppose that lim^ 00 ((yf) exists. Since d(g)/dL < [Eq. (|5.2|) 1 this implies 
lim^oo d{g)/dL = 0. Hence limL^oo^ 2 ) = by Eq. (]5^) . Since (g) 2 < (g 2 ) this im- 
plies that also limx,_ >00 (p) = 0. 



Where does the non-zero limit in Refs. [lCj and [TT| come from? The ground-breaking 



contribution of Zirnbauer was to use a "super- Fourier expansion" of the heat kernel W(Q, Q 1 ) 
in terms of eigenfunctions of the Laplacian in the space of the a model. This type of 
Fourier analysis is well understood for classical symmetric spaces. The development and 
application of the supersymmetric analogue for the a model enabled Zirnbauer, Mirlin, and 
Miiller-Groeling to compute non-perturbatively the first two moments of the conductance 
for any (3. The non-zero limiting value lirrii_ KX) (p) = 1/4 for (3 = 4 resulted from a "zero 
mode", a non-trivial eigenfunction of the Laplacian with zero eigenvalue. Since this zero 
mode does not decay as L — > oo, it led to the surprising conclusion of absence of localization 
in a wire with spin-orbit scattering in zero magnetic field.0 



An explicit expression for the zero-mode was not obtained in Refs. |10| and |Tl|, but only 
its contribution to the moments of the conductance was computed. By inspecting the initial 
condition ( |3.15| ) of the generating function for the a model we have been able to construct 
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a zero mode for /3 = 4, but only if we ignore the Kramers degeneracy of the transmission 
eigenvalues. This unphysical zero mode, given by 



^ ov ' 3 ' ; 2 + 2cosh^i ' v ; 

arises by taking the initial condition fl3.15|) with Ni = 1 and (3 = 4, but without the 
extra factor two in the exponent, required by Kramers degeneracy. This unphysical initial 
condition solves the evolution equation ( |3.11| ) for the ensemble average of the generating 
function and implies an L- independent average conductance (g) = 1/4. Although we can 
not prove that Eq. Q5.3| ) is Zirnbauer's zero mode, the coincidence with the limiting value 
lini^oo^) = 1/4, linif^oovarg = 1/16 is quite suggestive. 

The reason why we have to exclude the zero mode ( |5.3|) from the Fourier expansion of 
the heat kernel is that it is not single-valued on the a model space of supermatrices Q, 
although it is a well-defined function of 9. The parameterization ( p.l|) of Q is 27r-periodic 
in the angles 9± = 9 3 ± 6*4. We can then consider on the space of angles 6*3, 6*4 a parity 
operation P which consists of adding tt to both of these angles. This parity operation does 
not change Q, but it changes the zero mode ( |5.3j ). The Laplacian ( |2.4b[ ) commutes with P 



and the eigenfunctions have therefore either even or odd parity (eigenvalues +1 or —1 of P, 
respectively). The physical modes of the a model must have even parity, since only these 
functions are single-valued. For (3 = 4, it is the Kramers degeneracy which ensures that the 
initial condition fl3.15|) has even parity. 



This observation led us to check the parity of the eigenfunctions 4>v(Q) of the Laplacian 
in the super Fourier analysis of Refs. [l^ and [D] We consider the eigenvalue equation 



A^„(0!, 9 3 , 4 ) = -e(i/)&,(0i, 9 3 , 4 ) (5.4) 

for (3 = 4 in the limit 9\ — > 00 at fixed #3, 64. In this limit, the Laplace operator simplifies 
considerably 

a a, 9 9 d 1 d . _ d 1 d . _ d 
Wi Wi sm #3 afc's c/6'3 sin 4 0^4 a0 4 

From this expression one may identify the set of quantum numbers v — (A, 1 + 2ni, 1 + 2n2), 
where A is a real number and ni, ri2 are non-negative integers. The asymptotic behavior of 
the eigenfunctions 4> v {6\, 9 3 , 9 4) is given by 

Mdi , O3, 9 4 ) ~ exp \\ (1 + iA)0il (P ni (cos 9 3 )P n2 (cos 9 4 ) + P m (cos 9 3 )P m (cos 4 )) , (5.6) 



with the Legendre polynomials P n (x) and the eigenvalues 

e(A, 1 + 2n x , 1 + 2n 2 ) = - (^A 2 + (1 + 2m) 2 + (1 + 2n 2 ) 2 - l) . (5.7) 

The parity of this eigenfunction is just (— l) ni+ ™ 2 and we have to restrict ourselves to those 
Hi an d n 2 with rii + n 2 even. Applying this selection rule to the expressions for (g) and (g 2 ) 
of Refs. |TD| and |TT|, omitting the zero mode [and the subsidiary series with quantum numbers 
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v = ± 2) of Refs. |T^ and [11], for which the asymptotic behavior ( |5.6| ) is also valid], 

and multiplying the surviving terms with a factor of 2 to account for Kramers degeneracy, 
yields for (3 = 4 and in the limit N% = N 2 = N — > oo the expression 



(g n ) 



)5— n 



X 



Y, / dX A(A 2 + l)tanh(7rA/2)ZiZ 2 p„(A,Zi,Z 2 ) 

i 1 ,I 2 = l,3,5,..., ^° 
i 1 + i 2 =2 (mod 4) 

J| (-1 + ia\ + (Ti/i + cr 2 Z 2 ) _1 exp 



-(A a + + l)L/(2fl 



cr,cri,(T2=±l 



where n — 1,2 and 

Pi(A,/i,/ 2 ) 
p 2 (A,Zi,Z 2 ) 



A 2 + Z 2 + / 2 -l, 

i[2A 4 + Z 4 + Z 4 + 3A 2 (Z 2 + Z 2 ) - 2A 2 + Z 2 + Z 2 - 2]. 



(5.8a) 



(5.8b) 
(5.8c) 



Note that in our notations the dimensionless conductance g is by a factor 2 smaller than 
g in the notations of Ref. [□]. Comparison of Eq. ( |5.8|) with the f3 — 4 result of Ref. TT| , 
where the parity selection rule was not implemented, shows that the perturbation expansion 
around L/£ = is the same. (We checked this numerically up to order (L/£) 3 .) Outside 
the perturbative regime, the two expressions are completely different. Instead of a non-zero 
limit (g) = 1/4 for L/C, ^> 1, we find from Eq. ( |5.8| ) the exponential decay 



(5.9) 



(g) « f (2L/nO- 3/2 e- L ^. 



To test our result, we have compared it with a direct numerical simulation of the IWZ 
model by Mirlin and Muller-GroelingEll (with M = 100, iV = 25 and an average over 100 
different samples). The comparison is shown in Figs. ^ and |3|. It is clear that our Eq. 



( |5.8|) (solid curve) agrees quite well with the simulation, while the result of Ref. |TT] does not 
(dotted curve). 

Notice that this issue of the parity of the eigenfunctions does not occur for (3 — 1, 2, since 
there is only one compact angle (83) in those cases. The parity operation on the ^-matrices 
exists only for /3 — 4. For completeness we collect in Figs. |4] and [5] the results for (g) and 
varg for all three symmetry classes. The (3 = 1, 2 results are from Ref. 0, the (3 = 4 result 
is our Eq. (|5T 



VI. CONCLUSION 

We have established the exact mathematical equivalence of the two non-perturbative 
theoretical approaches to phase-coherent transport and localization in disordered wires: 
The Fokker-Planck equation of Dorokhov, Mello. Pereyra, and KumarEH] and the one- 
dimensional supersymmetric nonlinear a model.i@'i0ii The equivalence has the logical 
consequence that the absence of localization in the symplectic symmetry class, obtained 
by Zirnbauer by super- Fourier analysis of the a model, is not correct. By applying a selec- 
tion rule enforced by Kramers degeneracy to the eigenfunctions of Refs. [K] and [11], we have 
obtained modified expressions for (G) and var G, which decay exponentially as L — ► 00 and 
which agree well with existing numerical simulations.ill 
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Our equivalence proof has both conceptual and practical implications. The DMPK 
equation and the Id a model originated almost simultaneously in the early eighties, and 
at the same institute.01 Nevertheless, work on both approaches proceeded independently in 
the next decade. Knowing that, instead of two theories, there is only one, seems to us a 
considerable conceptual simplification of the field. It implies that the microscopic derivations 
and random-matrix models developed for the a model apply as well to the DMPK equation, 
and vice versa, (we see only the restriction, that the a model requires the thick-wire limit 
N — > oo, while the DMPK equation applies to any number of channels N.) Practically, 
each of the two approaches has its own advantages, and now that we know that they are 
equivalent, we can choose the approach which is best suited to our needs and skills. 
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APPENDIX A: TRANSPORT PROPERTIES DETERMINED BY THE 

GENERATING FUNCTION 

We list the transport properties of interest that can be generated from F(9, A), following 
RejaeiH Let us consider the function 



/det(l + z 2 tH) 
f{ZuZ2) -\det(l + Zl tH) 



(Al) 



which equals (F(6,X)) at z 2 = — sxr 2 (^9 s ), Z\ = sinh 2 (|6' 1 ), and 9 2 = 9^ = 0. We write Eq. 
~TT) in the form 



f{zi,z 2 ) 



exp 



Tr ln(l + z 2 tH) - Tr m(l + z x tH) 



(A2) 



The standard expansion™! with respect to small z\ and z 2 yields the first two moments of 
the dimensionless conductance g = (l/d)TrtH (with d = 1 + S^p), 



{9) = \-K~f(z 1} z 2 ) 

d OZ 2 zx=0=z 2 

1 9 f< \ 



(g 2 ) 



1 d d 

d? dz\ dz^ 



21=0=22 



(A3) 
(A4) 
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We may also consideJllE derivatives of f(zi,z 2 ) at z\ = z%- This may require the 
analytic continuation of 9\, 63 to complex values if Z\ < 0, z 2 > or z 2 < —1. Therefore, we 
introduce the function f{z\) as 



a 



f(zi, z 2 ] 



dz 2 

(Tr [(1 + z 1 tH)- 1 th 

00 

n=0 

^Wl) - (Tr [(l + Sit+t)- 1 !)). 



The average density of transmission eigenvalues now follows from: 



p(T) = (Tr6(T-th)) 

= --^Im/(-(T + ^0 + )- 1 



(A5) 



(A6) 



The application of Eq. ( |A6[ ) requires the analytical continuation of both variables z\ and z 2 
to values < — 1. 



APPENDIX B: THE GENERATING FUNCTION IN TERMS OF THE 

TRANSMISSION MATRIX 



In this appendix, we show that Eq. ( |3.5| ) for the generating function in the IWZ model 
equals Eq. ( |2.3| ). We first consider the two cases /3 — 1, 4 of time reversal symmetry, when 
7i = Hl 8 in Eq. ( |3.5a| ). The necessary modifications for (3 = 2 are described at the end. 

We make use of the folding identity 

Sdet fj = Sdet(l„ - AB) = Sdet(l m - BA). (Bl) 

We abbreviate G± = (E -H ±mWW^)~ 1 . Taking out the factor (E -H + iirWW^A) (with 
unit superdeterminant), we may rewrite Eq. (|3.5a|) as 

F = SdeH ( 1 + ( G + £ j zW^Q - A) J 



^G+Wi 
WlG-W x 



Sdet-Ml + z 7 r(^^^ T " , ](g-A)], (B2) 



where we have applied Eq. (|B1|) on 5 = W^. We now insert the reflection matrix 
1 - 2mWlG + W 1 [see Eqs. and Q] into Eq. Q and obtain 



F = SdeH Q(l + AQ) + i f J r ° t ) (1 - AQ)) . (B3) 
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Now we use the parameterization (|2.1 ) for Q. Notice that Eq. ( |B3| ) does not depend on the 
angular part of Q [the matrices u, v in fl2.1b|) 1. Hence we may choose Q as 



Q = T- 1 AT, T = exp(±0), 9= f ° ^ 



(B4) 



which leads to 

|(1 + AQ) = cosh(ie) T, |(1 - AQ) = - sinh(|e) T. 
Inserting this in Eq. ( |B3| ) and taking out the factor T (with unit superdeterminant), we get 

F = Sdet _ 1/2 / cosh(I^) -rsinh(§0)\ 
V — rt sinh(|^) cosh(i^) J 

= Sdet' 1/2 (cosh 2 (±fl) - sinh 2 (i0)rV) 

= Sdet~ 1/2 (l + sinh 2 (±fl)^) , (B5) 

where we have again used flBl|) and the relation r*r = 1 — vt imposed by unitarity of the 
scattering matrix. 

The matrix vt has eigenvalues T n = (1 + A„) _1 (n = 1, . . . , Ni), which are twofold 
degenerate for (3 = 4, hence 

F = J] Sdet- d/2 (l + sinh 2 (±#) T n ) 

n=l 

= [] Sdet" d/2 (A n + cosh 2 (±#)) . (B6) 

n=l 

Here d — 1 + Eq. ( |B6| ) establishes the connection between the two expressions ( |2.3|) 
and (|3.5| ) for the generating function. 

The above calculation assumes Tt = His, hence (3 = 1 or (3 = 4. In the case /3 = 2, 
one has instead 7i = Re(i7)l 8 + i Im(if )r 3 . Instead of two subblocks with r and [see Eq. 
(|B3|)1, one now needs four subblocks with r, r T , r* , and r*. Repeating the calculation one 
finds that the final result (IB 61) remains the same. 



APPENDIX C: IDENTITY OF LAPLACIANS 

The goal of this appendix is to prove Eq. (|2.6|) . Hereto we first analyze the structure of 
the l.h.s. of Eq. (|2.6| ) in more detail. 

The derivatives of F(9, A) with respect to 9j are calculated using 

~ IS 70m ^0- ) p v> A) ' (C1) 

~ IS mx n) 00] ) F ^ A) + 

df{§,X n )df(B,X r 



^ ^ x x / S ^ A) ' (C2) 

K n^mf(0,\ n )f(0,\ m ) °Vj <™j 
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Inspection of Eqs. ( |2.4b| ), (|C1|) and (|C2|) shows that AqF(9, A) has two contributions, one 
involving a single summation over the channel indices n, and another one involving a double 
summation over channel indices n ^ m, 



N 



N 



A § F(9, A) = [ £ ^(0, A n ) ] F(6, A) + ( £ <? 2 (0, A n , A m ) ] A). (C3) 



n=l 



Using the definition ( 2.3b|) of f(9, A) one may straightforwardly calculate the functions gi 
and g 2 . The expressions are rather lengthy and will not be given here. 

The r.h.s. of Eq. ( |2.6|) contains the differential operator D~ x , which is given by Eq. (|2.5b[) . 
Simple algebra yields 



N / q2 a 

D x = <r[\ n (i + \ n )— + (i+2\ r ; ' 



n=l 



d\ n . 



nJ=m \ n 



+ 



Am(l A m 



d 



OK 



As a consequence, D^F(9, A) has again the structure of Eq. 
by 



<?i(M) 
g 2 (9,Xi,X 2 ) 



1 



/(M) 
Ai — A 2 



A(l + A) 



d 2 /(M) 



+ (1 + 2A) 



(C4) 

with gi and g 2 now given 



dA 2 v ' dX 
X 1 (l + X 1 )df(9,X 1 ) X 2 {l+X 2 )df{9,X 2 ) 



/(Mi 



<9Ai 



/(Ms 



<9A 2 



(C5a) 
(C5b) 



Comparison of Eqs. (|C3| ) and ( p5|) shows that the two definitions of the functions g\ and g 2 
are identical. This completes the proof of Eq. (|2.6|) . 



APPENDIX D: EXTENSION TO HIGHER DIMENSIONAL SUPERMATRICES 



The argumentation presented in Sees. |TTT| and [TV] can be generalized to a models with Q 
matrices of arbitrary dimension 8p with p > 1. This generalized equivalence proof applies 
to the p-point functions p p (Ti, . . . , T p ) instead to the limited number of statistical quantities 
that can be generated by the "small" a model with p = 1 (compare appendix |A]). Here, we 
briefly present the modifications with respect to the p = 1 case. The modifications concern 
the parameterization ([2.1|) and the generating function (|2.3b|) . 

The main technical difficulty in such a generalization is due to the radial part of the 
Laplace operator. The procedure to calculate it on conventional symmetric spaces is 
standard" and is carried over to the supersymmetric a models as described in appendix 
B of Ref. 0. It is now more convenient to use a slightly modified form of the parameteriza- 



tion of the Q-matrices, where 9 in Eq. ( |2.1a| ) is fully diagonal (rather than block-diagonal): 




(£) 



nm -^n^nmi 



yJnm, l<n,m<2p. 
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The symmetry restrictions are [cf. Eq. fl2.1d )l 

Vi+p 



Vi 

X i 
J. ; 



•^i+pi Vi 
X{-\-p 



if/? 

Vi+p if P 
if/? 



;d2) 



for i = In the case p — 1, we have the relations X\ = 9i + 9 2 , x 2 = Q\ — 9 2 

Ui — &s + $4, and 2/2 = 6*3 — #4 between these parameters and the 0j used in (|2.1d ). 



We can directly apply the results of appendix B in Ref. [11] which are given in terms of 

the so-called roots a(Q) [with given as in QB4j) 1 . The roots are linear functions of G which 

are the eigenvalues of the linear mapping ad(0)(X a ) := [0,X a ] = a(Q)X a , defined on a 

certain super-Lie algebra.0 The eigenvectors X a of the mapping are the root vectors, which 

do not depend on 6. The radial integration measure J{9) in Eq. ( |2.4b| ) can be expressed 
„|o| 



J{6) = llsinh m "[±a(e)], (D3) 

«>o 

where the integer m a is the multiplicity of the root a (the dimension of the root space). 
Both positive and negative values of m a can occur. The factor \ is due to the difference 



between the normalization ( p. la] ) of 6 and the one used in Ref. |TT|. In appendix A of Ref. 
|LT| ) explicit formulas for the roots as well as for the root vectors are given for the case (3 = 1, 
p = 1. 

We have calculated the roots and the root vectors for all f3 and arbitrary dimension 
8p of the Q matrices. For simplicity, we only present the results for the roots and their 
multiplicities. Let us denote with p x (p y ) the number of independent Xi (yi) parameters, i.e.: 
p x = 2p,p,p and p y = p,p,2p for (3 = 1,2,4, respectively. Note that (3p x = 2p y . We find 8 
different types of (positive) roots: 



a(0) 




(i <j < i<p x ), 




m a 


= P, 


a(6) 


= Xj - iyi 


(i < j < p x , i<i< 


Py), 


m a 


= -2, 


a(6) 


= i(yj - vi) 


(1 <j < l<Py), 




m a 


= 4//3, 


a(6) 


= 2xj 


(l<j<Px), 




m a 


= 13-1, 


a(6) 


= 2i Vj 


(l<J<Py), 




m a 


= 4/P-l 


a(9) 


= Xj + Xi 


(1< j < l<Px), 




m a 


= 13, 


a(6) 


= xj + iyi 


(l<j<Px, l<l< 


Py), 


m a 


= -2, 


a(Q) 


= i (yj + yi) 


(1< j < l<Py), 




m a 


= 41(3. 



(D4) 



The radial part of the Laplacian takes the form 

d d 



P.v 



dxj 



dxj 



3 Pv 



d_ 

dyj 



d_ 



(D5) 



The expressions ( ^.3a| ) and ( ^.5a| ) for the generating function now remain valid with the 
modified 9 of Eq. ( pi| ), and with Eq. ( |2.3b| ) replaced by 



Py 



P.r 



f(0, A) = + 2A + cosfoi)]]!! 1 + 2A + cosh(; 



-P/2 



(D6) 



i=l 



i=l 



18 



It is convenient to use the variables Uj = sinh 2 (|xj) and = — sin 2 (|?/j) in terms of which 
the Laplacian has the form 



Px C) O O Py O O 

A» = £ J' 1 t%«,a + «,) /,„ (j - fg^g^-yd + «,) (DTa) 

Pa: Py 

j= n n («*-^) v/, nii(«i-*i)~ a ( D7b ) 
><nNi+^)^ 1 iiHi+^)) 2/H . 

i=l i=l 

The generating function F(9, A) is given by 

N / Py Px \ 

^ A) = J] IK 1 + A « + + A " + «i) _/,/2 • (D8) 

n=l \i=l i=l / 

We have verified that the identity of Laplacians [Eq. [2.6[) 1 remains true for the modified 
expressions ( P7a| ) and (P8|) . The calculations goes in a similar way as shown in App. for 
p = 1. Now, we have to keep track of 7 different types of contributions with double and 
triple sums over functions of A n , Ui, Vj. 

In App. [A] we have shown that the average density of transmission eigenvalues p(T) 
can be obtained from the generating function ( |2.3a|) . Using the corresponding function for 
the higher-dimensional a model considered here, it is straightforward to get the p-point 
correlation functions p p (Ti, . . . , T p ). 
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FIGURES 



FIG. 1. Schematic drawing of the disordered wire and the leads according to the IWZ model 
description. The left lead (lead 1) contains N\, the right lead (lead 2) N% propagating channels. 
The number of propagating channels in the disordered wire is N. In the IWZ model, the disordered 
wire is divided into K segments, each having a random Hamiltonian drawn drom the Gaussian 
ensemble. The derivation of the Id a model from the IWZ model assumes N 3> 1, but allows for 
finite N\ and A^. 

FIG. 2. The average conductance (g) multiplied by 4L/£ = L/Nl for the symplectic symmetry 
class as a function of 4L/£ for N 3> 1. Shown are our result Q5,8| ) (solid), the numerical simulation 
of Ref. H| (M = 100, N = 25) (dashed), and the result of Ref. (dotted). 

FIG. 3. As in Fig. || for the variance var g of the conductance. 
FIG. 4. The average conductance {g) multiplied by f3L/£ = L/N£ for the three symmetry 



classes as a function of /3L/£ for 2V S> 1. The curves for (3 = 1,2 are taken from Refs. 10 and 11 
and the curve for j3 = 4 is calculated from Eq. (|5.8|). Notice that ^ = j3Ni is proportional to 0, so 
that the scaling of the axes is /3-independent. 

FIG. 5. Same as Fig. H, for the variance var g of the conductance 
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